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We consider the possibility of time variations of the solar neutrino flux due to 
. the radial motion of the Earth and neutrino interference effects. We calculate 

the time variations of the detected neutrino flux and the extent to which they 
are suppressed by Coulomb collisions of the neutrino emitting nuclei. 
CO ■ We treat the process of neutrino emission under the second quantization 

■ method. A neutrino emitting nucleus suffers random Coulomb accelerations 

in different directions. Thus, the positions and the velocities of the nucleus 
are different at different emission times. We properly calculate these emission 
times and emitted neutrino energies by making use of the method of sta- 
tionary phase, which leads to the causality and Doppler shift equations. We 
average our results over the Holtsmark distribution of nuclear accelerations. 
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O { To properly treat the collisions, it is necessary to simultaneously include in 

our analysis all other significant physical decoherence effects: the energy av- 
eraging and the averaging over the position of neutrino emission. 

We find that the collisional decoherence is not important (in comparison 
with other decoherence effects) in the case of neutrino line, but it may be 
important in the case of the continuous neutrino spectrum (or for high energy 
resolution explorations of the internal structure of the line). We find that 
the collisional decoherence averages out the time variations for the neutrino 
masses m\ — m\ > f -7 eV 2 S^eV- 

A simple and clear physical picture of the time dependent solar neutrino 
problem is presented and qualitative coherence criteria are discussed. Exact 
results for the detected neutrino flux and its time variations are obtained 
for both the case of a solar neutrino line, and the case of the continuous 
neutrino spectrum with a Gaussian shape of the energy response function of 
the neutrino detector. We give accurate constraints on the vacuum mixing 
angle and the neutrino masses required for flux time variations to not be 
suppressed. 

Pac(s): 26.65.+t, 14.60.Pq, 96.60.Jw 
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I. INTRODUCTION 



To explain the discrepancy between the observed flux of solar neutrinos and the predicted 
emission rate it has been postulated that the neutrinos have two or three different mass 
states, and that these mass states are different from the flavor states This extension 

of the standard electroweak theory leads to oscillations between the electron neutrino and 
the different neutrino species. A key role in the solution of the solar neutrino problem is 
played by the MSW theory |8|||. In this theory there is a spherical resonance region at the 
radius r = r res . We denote the region inside the resonance radius, r < r res , as region / and 
the region outside the resonance radius, r > r res , as region II (see Fig g). Hereafter, we 
suppose the solar core to be inside region I (i. e. inside the resonance radius r = r res ). In 
this case all neutrinos are emitted in region I and they cross the resonance region once to 
enter region II on their way to the Earth (as shown in Fig |2|). The other situation, when 
the resonance region in inside the solar core, so that neutrinos emitted in region II cross 
the resonance twice or not at all, presents complications in notations rather than in physics, 
and it is discussed in Appendix A. 

In our paper, we consider the case of two types of neutrino flavors and correspondingly 
of two types of neutrino mass states. In this case for each neutrino energy there are two 
orthogonal WKB solutions of Dirac equation out of the resonance region, which we denote 
as + and — neutrino states. These neutrino states exist in regions I and II, where the WKB 
approximation is valid, and they coincide with neutrino pure mass states in the vacuum, but 
they are not given for the resonance region, where the WKB approximation breaks down. 

In the adiabatic approximation the + neutrino states are connected to each other as well 
as the — states are connected. So, we say that there are two neutrino paths: (+, +) and 
(— , — ) [for example, the (+, +) path means that the neutrino is emitted in region I as the + 
state and is converted, while crossing the resonance, into the + state in region II] . However, if 
the adiabatic approximation breaks down, then the + and — states jump one into another at 
the resonance radius. As a result, in this case there are four neutrino paths: (+, +), (+, — ), 
(— , — ) and (— , +) (see Fig |]). Due to electron scattering and nonzero neutrino masses 
these neutrino paths gain different quantum mechanical phases during the propagation to 
the Earth and hence interfere with each other. 

There are two neutrino eigenf unctions for each neutrino energy, which are solutions of 
Dirac equation in the entire region of space (including the resonance region, where the 
WKB approximation is not valid). It is convenient to choose them so that they coincide 
with the + and — states in region I (where neutrinos are emitted). In region II (after the 
resonance crossing) the eigenf unctions are sums of the + and — states, whose coefficients 
are determined by resonance jump formulas ( |3.6| )-( pT77| ) [see Fig 0j. For this choice the entire 
eigenfunctions are orthogonal. Each neutrino propagating to the Earth with a given energy 
can be represented as a linear combination of two eigenfunctions or as a linear combination 
of four neutrino paths (with only two independent coefficients and all others related to the 
jump formulas). For example, in the adiabatic approximation there is no jump, so, there 
are only two paths, the (+, +) and (— , — ), and the two eigenfunctions correspond to these 
two paths (out of the resonance region). 

The MSW effect results in a reduction of the average number of detected electron neu- 



trinos [10] and also produces time variations of the neutrino flux. The time variations of 
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the detected neutrino flux are due to interference between different neutrino paths and are 
sensitive to the neutrino masses and the mixing angle. There are several physical effects 
that can destroy the coherence between neutrino paths and limit the observability of flux 
time variations. The most significant effect is an averaging of the observed neutrino flux over 
the effective detected energy band a. This averaging wipes out a cross term between two 
neutrino paths if a > at == h/r, where r is the difference of propagation times of these paths. 
For a neutrino line, a is simply the energy width of the line, a ~ T (thermal broadening 
effect). In the case of a detection of the continuous neutrino energy spectrum the energy 
response function of our detection device must be sufficiently narrow if one wants to observe 
flux time variations. 

A second effect, which can destroy coherence, and which is equally important, is the 
averaging of the detected neutrino flux over the region of neutrino emission inside the sun. 
This effect depends on the electron density in a region of the solar core and it is important 
even for the case of very small neutrino masses (see Sec. [V Aj) . As we will see below, the 
averaging over the position of neutrino emission greatly simplifies all calculations because 
only two cross terms between four neutrino paths survive. (Actually one of these cross terms 
is much larger than another, so there is only one potentially important cross term! This 
cross term varies with the motion of the Earth and therefore it gives time variations of the 
detected neutrino flux.) 

A third effect that can potentially average out time variations of the neutrino flux is colli- 
sional decoherence. Electric forces acting on the emitting nucleus from its neighbors lead to a 
negligible Stark effect on a nuclear scale but they also lead to the acceleration of the nucleus. 
The collisional decoherence was first discussed by Nussinov in 1976 [IT], who used a coherence 
criterion r < h/AS with AS approximately equal to the electric potential due to neighbor- 
ing ions. For the 0.862 MeV 7 Be neutrino line Nussinov got r < 3 x 10~ 17 sec, although the 
theoretical basis is uncertain. In 1985 Krauss and Wilczek |12| suggested the plasma full 
collision time as an approximate value of the coherence time. They got r < 10 _15 sec. As we 
will see, this is a considerable overestimate for the coherence time. More recently, Loeb [TJ| 
considered a collisional Doppler spectral broadening (|A£|) = £{\AV\)/2c. Loeb obtained 
a criterion r < (2X/a) 1 ^ 2 , where a = (l/W^) -1 was the harmonically averaged nuclear ac- 
celeration and v co \ was the collision frequency. Loeb's coherence time for the 0.862 MeV 7 Be 
neutrino line was 5 x 10~ 17 sec. 

Because of the uncertainty in these estimates of the criterion for the collisional deco- 
herence, we felt a more detailed analysis was desirable. In this paper we carry out such an 
analysis and derive the correct collisional decoherence criterion. It turns out that Nussinov's 
and Loeb's estimates are close to this criterion if the acceleration they used is replaced by 
the Holtsmark acceleration. 

In most of the previous papers the behavior of the neutrinos is treated as a one- 
dimensional problem, where a constant neutrino energy is assumed and the evolution of 
the neutrino in space (or equivalently in time) is followed as it propagates to the Earth. 
The neutrino is assumed to be emitted into a single state in the high density solar core. It 
propagates to the resonant region where it nonadiabatically splits into two states. These 
two propagate with different momenta to the Earth, and at the detector on the Earth these 
two states can interfere according to the accumulated phase difference during this latter 
propagation. 
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This procedure is adequate for treating the decoherence due to the energy averaging. 
However, it is not adequate to properly treat the collisional decoherence of the neutrino 
flux time variations, since the neutrino emitting nucleus moves erratically in time (suffering 
random accelerations in different directions). Thus, the positions and the velocities of the 
nucleus are different at different neutrino emission times. These emission times have to be 
properly defined and evaluated. To do this, both space and time must be introduced into 
the problem. 

In our calculations we follow the prescription of Fermi. We assume two definite spatial 
eigenfunctions for each neutrino energy (because we assume two neutrino flavors). We then 
expand the neutrino wavefunction in terms of these eigenfunctions. The components of this 
expansion are functions of time alone, as the neutrino is emitted by the moving nucleus. The 
components themselves have probability distributions, and after second quantization these 
distributions are, in turn, expanded into two oscillator states corresponding to occupation 
or non-occupation of the neutrino eigenfunctions. Finally, the differential equations in t and 
x for the time evolution (during emission) of the coefficients a of the latter expansion (i. e. of 
the oscillator occupation amplitudes) are derived and integrated. 

For a neutrino line the decaying nucleus emits a neutrino with the fixed energy in the 
frame co-moving with the nucleus. By integrating the differential equations one finds that 
each of oscillator occupation amplitudes a is an integral over a time parameter if correspond- 
ing to "the time of emission" . Then we sum over excited oscillator states by integration over 
a spread of neutrino energies £ in the laboratory frame and summing over ± indices (the 
summing over ± is actually the summing over two eigenfunctions for each energy, which 
coincide with the + and — adiabatic states at the position of emission). Thus, the emitted 
neutrino wavefunction is a double integral over t' and £ and, of course, a function of t and 
x. Because of the rapidly changing phase the integrands in this double integral vary rapidly 
with t' and £. For a given time and position of detection t and the main contribution 
to the integral comes from isolated values for t' and £ where the phase is stationary. As 
one would expect, these values coincide with the solutions of the causality and Doppler 
shift equations for the emission time t' and the emitted neutrino energy £. The method of 
stationary phase selects these solutions and also gives the neutrino phase. 

Let us, for simplicity, consider that the neutrino is emitted into only one of two neutrino 
eigenfunctions, for example, into the one that coincides with the — state in region I (as 
in case of high electron density in the solar core). Then, in region II this eigenfunction is 
given by the sum of — and + states due to nonadiabatic resonance crossing. In other words, 
the entire eigenfunction is the mixture of two neutrino paths (— , — ) and (— , +) out of the 
spherical resonance region (i. e. in regions I and II). These two neutrino paths have the 
same momenta fc_ in region I but different momenta and k+ in region II [k- < k+, see 
Eqs (p.lO| ) and (|2.12Q 1. Therefore, for a given time and position of the detection, one finds 
two solutions for t' and £ from the method of stationary phase, which correspond to the two 
different propagation times of the two neutrino paths. We can say that the eigenfunction 
is emitted twice: first the (— , — ) path contributes, then the (— , +) path contributes. The 
two different values of the emission time t' are determined as the two retarded times. The 
corresponding two different values of the neutrino energy £ are found to be the two Doppler- 
shifted energies, since the emitting nucleus has two different velocities at the two emission 
times due to its acceleration. In terms of these parameters, t' and £, the relative phase of 
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the two neutrino paths (i. e. of the two solutions for the emission of the eigenfunction) can 
be found and the amount of interference between them can be determined. Then we have to 
average over the distribution of nuclear accelerations, since we detect many emitting nuclei 
and integrate over a large interval of detection time. Only in this way can we correctly 
determine the decoherence produced by collisions. The resulting decoherence depends on 
neutrino masses and the mixing angle. 

For the continuous neutrino spectrum one needs to appreciate that an electron is emitted 
as well as a neutrino. To carry out the problem in this case one, first, fixes the state into 
which the electron is emitted and then treats the neutrino emission as in the case of neutrino 
line emission. However, at the end a sum over electron states must be taken. This introduces 
some further complications. For example, the two neutrino paths have different energies, so 
if the energy band width of the detector is too narrow, the detector will be affected only 
by that one of these neutrino paths which has its energy inside the detector band, and the 
other path will be outside this energy band. Thus, no interference can occur. In the case of 
the continuous neutrino spectrum detection both collisional decoherence and averaging over 
the neutrino energy must be treated simultaneously. 

In this paper we present a complete self-consistent analysis (making use of the method 
of stationary phase) of the time variations of the detected neutrino flux due to interference 
effects. We consider a general nonadiabatic theory for two neutrino flavors. We neglect 
spherical geometry of the problem and use jump formulas derived under the assumption 
that the solar electron density varies exponentially with the radius throughout the resonance 
region. We make use of the last assumption only for our final numerical results. Another 
assumption that we make is that of constant nuclear acceleration, i.e. the acceleration of 
the emitting nucleus is considered constant during the delay time r between propagations of 
different neutrino paths. However, the nuclear acceleration has random direction. We use a 



Holtsmark distribution for the absolute value of the instantaneous nuclear acceleration |L4 
The assumption of constant nuclear acceleration turns out to be sufficient for our purposes, 
i.e. it is valid for all relevant cases in which time variations of the detected neutrino flux are 
significant. 

In the second section we use the WKB approximation to solve the Dirac equation in the 
mass state representation and to derive formulas for the + and — neutrino states. Then, 
in Sec. [TTI| we discuss the nonadiabatic jump at the resonance radius and we write down 
the two neutrino eigenfunctions. In Sec. |IV A| we study the processes of neutrino emission 
and detection under the second quantization theory. Using the method of stationary phase, 
we obtain the neutrino wave function at the Earth and the causality and Doppler shift 
equations in Sec. [IV B|. These equations are solved and the quantum mechanical phase 



differences between neutrino paths are calculated in Sec. |IVC| making use of an expansion 
in small neutrino masses. Since the final exact calculations are quite complicated, we give 
a qualitative analysis of the problem in Sec. [V] (which is devoted to the averaging over the 
effective detected energy band and to the collisional decoherence) and in Sec. |V A] (which 
discusses the averaging over the region of neutrino emission inside the sun). The coherence 
criteria and some of the numerical results are summarized in Fig. [| In Sec. [V| we derive 
the exact numerical results for time variations of the detected neutrino flux. Section |V1 A 
with Figs. |4], [5] presents the exact analytical calculations and numerical results for the case 
of a solar neutrino line. Section |V1 Ij| with Figs. |5|-^ presents the same for the continuous 
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neutrino spectrum. For the continuous neutrino spectrum we assume a Gaussian energy 
response function for our detection device. In Sec. VU we give the conclusions. Finally, in 
Appendix A we discuss a generalization of the concept of neutrino paths, and in Appendix B 
we give the formal mathematical proof of the averaging out of most of the cross terms between 
different neutrino paths due to summing over the region of neutrino emission. Analytical 
evaluations of some integrals are given in Appendices C and D. 



II. THE WKB STATES 



Let there be two mass states of the neutrino with masses mi and m 2 (we assume that 
m 2 > mi). Then one writes 



\b> e > = cos^lz/i > + sin^l^ >, 
\v x > = — sin 6\ui > + cos^lz^ >, 



(2.1) 
(2.2) 



where u e refers to the electron flavor and u x to the other flavor, possibly the muon flavor. 9 
is the vacuum mixing angle (we assume that < 9 < 7r/4). Further, the electron neutrino 
flavor interacts with the electron density n e of the sun more strongly than the | v x > flavor, 
so that the \v\ > and \u 2 > states are coupled. (We neglect neutral-current interaction 
which is the same for both flavors.) Each of the mass states for a neutrino propagating in 
the x direction is a two-component state consisting of the first and fourth components of 
the Dirac wavefunction. In an infinite homogeneous medium, with electron density n e , the 
Dirac equation, in the mass state representation, for the four components of the neutrino 
wavefunction, two for each mass state, is |jl0l,3ffl|l5H 



d 













j' x Vx + 1 m 1 







1 g cos 2 9 1 g cos 9 sin 9 




U 





j' x p x + 1 m 2 _ 


+ 


1 g cos 9 sin 9 1 g sin 2 9 


) 


V 



(2.3) 



== \f2G?n e describes neutrino interaction 
1. In vacuum, when g = 0, U is the 



where the matrix composed of the elements with g 
with electrons, p x = —id/dx, and we take % — 1, c 
two-component wavefunction corresponding to the mi state and V is the two-component 
wavefunction corresponding to the m 2 state. The two-component matrices 



It 



1 
-1 



and 



7, 



1 
-1 



(2.4) 



are the first and fourth components of the usual Dirac matrices -y t and ~y x , while 1 is just 
the two by two unit matrix JTI|. We treat the neutrino wavefunction as one-dimensional in 
space. 

Now, g depends on x since n e is nonzero and varies inside the sun, and g is zero be- 
tween the sun and the Earth. (We replace the radial coordinate by x and ignore spherical 
effects.) Consequently, the solution of the Dirac equation is space dependent. We write 
two-component wavefunctions XJ(x,t) and V(x,t) as those with the same energy £ 



U(x,t) 



ui(x, £) 
u 2 (x, £) 



exp(—i£t), V(x, t) 



vi(x,£) 
v 2 (x,£) 



exp (—i£t). 



(2.5) 
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In the limit that mi and ni2 are small compared to S, the solution of the Dirac equation, 
Eq. ( |2.3|) , with U and V given by ( |2.5|) is 



U 2 



U(x,£) 



1+mx/S 
1 



and 



V-2 



V(x,£) 



1 + m 2 /£ 
1 



(2.6) 



where the two spatial scalar functions U and V are determined at energy £ by 

0. (2.7) 









V 



p x — £ + ml /2£ + g cos 2 9 gcos9sm9 

g cos 9 sin 9 p x — £ + m\j2£ + g sin 2 9 

Moreover, we can regard g as a slowly varying space function over most of the distance 
between the source and the Earth with the exception of the spherical resonance region in 
the sun, where g ~ A = f (m 2 , — m\)/2£. Therefore, we take the space dependence of U and 
V in the WKB approximation [TH] as 



U M ^ M ^J kdx ), V^)^)^/^), (2,) 

where X-^ is the position of neutrino emission, and u(x, £), v(x, £) are slowly varying scalar 
functions of x. As a result, equation (|2.7|) becomes 



k — £ + m\j1£ + g cos 2 9 g cos 9 sin 9 

g cos 9 sin 9 k — £ + m\/2£ + g sin 2 # 



0. 



(2.9) 



The space dependence occurs in g(x). 

For every £ there are two solutions of this normal mode equation. Let 



k = £ 



(2.10) 



The determinant of the matrix of normal mode Eq. ( |2.9| ) is zero, so k' satisfies the equation 

(2.11) 



k' 2 + k'g(x) - ^A 2 + ^A • g(x) cos 29 = 0. 



Hence, the two values for the local momenta k'± are given by 



k'± = —~g ± 9 2 - 2Ag COS26 1 + A 2 = --g ± - \[{g - A cos 29) 2 + A 2 sin 2 29. (2.12) 
The corresponding solutions for the WKB region are given by the + and — neutrino states, 



\il>±> 



u± 
v± 



x 

exp ( — i£t + i J k±dx 



g sin 9 cos 9 
-k' ± + A/2 -g cos 2 9 



X 

exp ( — i£t + i J k±dx 



(2.13) 



or alternatively as, 
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u± 
v± 



exp 



i£t + i k±dx) = ±N± 



k' ± 



+ A/2 + #sin 2 # 
—g sin 9 cos 9 



exp 



k+ dx 



(2.14) 



where iV± and N± are positive normalization constants. Equations (|2.10 ) and (|2.12|) — ( 2.14|) 
give us two WKB solutions of the wave equation HJ4J]. Note that these + and — solutions 
(states) are orthogonal at each point. 

At the Earth, X E , when g = 0, we have from Eqs. (|2.12|) , (|2.10|) 



k' ± (X E ) = ±A/2, 



and 



k+{X E ) = £-m\/2£, k_(X E ) = £ - m 2 2 /2£. 
In this case equations ( |2.14j ), ( |2.13| ) give us 



~u + (X E )~ 




' 1 " 




'u_(X E )~ 




"o" 


v+{X E )_ 







> 


v-(X E )_ 




1 



(2.15) 



(2.16) 



(2.17) 



In other words, at the Earth 





= \ut> 


= cos 9 1 v e > 


— sin 9 


\ v - > 


= \u 2 > 


= sin 9 1 v e > 


+ cos9\ 



(2.18) 
(2.19) 

On the other hand, if we assume that at the position of emission, X^, g ^> A, then from 
the same Eqs. (gig) , (gig) and ( gig ) we have 



-A cos 20, fel(^ 



-p + -Acos26>, 
2 



and 



|i/+> 



(A/5-)sin6'cos6'|z/ e > - [l - (A 2 /2g 2 ) sin 2 9 cos 2 6»j |z/ x >, 
1 - (A 2 /2sf 2 )sin 2 6 l cos 2 6'l |z/ e > + (A/p) sin6» cos0\v x > 



(2.20) 



(2.21) 



In the adiabatic approximation the WKB regions are connected continuously across the 
resonance region, so the + states are connected and the — states are connected. As a 
result, for the adiabatic theory, in the case of g(X-$) ^> A, the emission of the neutrino is 
primarily into the — state, which at the Earth makes a small contribution of sin 9 (eq. |2.19|) 
to the amplitude for an electron neutrino detection (if the mixing angle is small), while 
the contribution into the + state, ~ (A/g) sin cos 0, is small, but this state makes a finite 
contribution, cos 6* (eq. |2.18|) to the amplitude for the neutrino detection. To clarify the 
statement above we plot the dependence of neutrino relativistic mass m 2 el = £ 2 — k 2 on the 
electron density for both states in Fig. |l] . The right side of the figure corresponds to high 
electron density in the sun interior, the left side to low electron density at the Earth. 
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III. NONADIABATIC EFFECTS AND NEUTRINO EIGENFUNCTIONS 



The solution of the Dirac equation, the + and — neutrino states, were found in the 
previous section under the assumption that the WKB approximation is accurate. However, 
for the resonance region, 



where 



\g-g VC s\ S A, 



g ves = A cos 29, 



(3.1) 



(3.2) 



the WKB approximation breaks down because k + and k_ are very close at the resonance 
point, 



k + — /c_ = k', — k'_ = A sin 29 for g = g r 



(3.3) 



[see Eq. (|2.12|) 1. As a result, the + state can jump into the — state and vice versa []lO|fl- 
One has to find the probability of the jump to connect WKB solutions for the neutrino 
wavefunction in the regions I and II. The jump probability 5 2 can be derived under the 
assumption that the electron density n e varies linearly throughout the resonance region |IO|fl 



5 2 



exp 



vr sin 2 29 
-A 

2 cos 29 



n,. 



\dn e /dr\ 



(3.4) 



or under the assumption that n e varies exponentially throughout the resonance region ||r7| . [18 
exp (2tt r A cos 2 9 



1 



exp (27rr A) — 1 



exp 


10 3 cos 2 9 £ M l Y (mj - m\] 


-6~ 


- 1 


exp 


103 ^Mev ( m 2 -ml)_ 6 _ 




1 



(3.5) 



Here r 



6.6 x 10 9 cm is the exponential density scale 

,2 _ z' 2 



[\. e. n e [r) 
-l 



n e (0)exp[r/r o ]), 

and we use notations (m 2 , — m'()_ 6 = (ml — m 2 )/10~ 6 eV 2 and ^Mev = ^/MeV. In our 
numerical calculations we use formula (|3.5| ) for the probability of the jump, because it is a 
good approximation for the whole range of neutrino parameters in the relevant case when 
the resonance electron density is much smaller than the density inside the solar core (see 
discussion and results below). 

Note, that for very small neutrino masses, ml — m\ 10 _9 eV 2 £]yieV; the jump prob- 
ability (|3.5| ) has the limit 5 2 — > cos 2 9. It is so-called "vacuum oscillations" case, when the 
neutrino is an oscillating electron neutrino in the vacuum. This extreme non-adiabatic limit 
is included into our calculations and results. However, we do not use the term "vacuum 
oscillations" in our paper, because we believe it could be misleading, and the solar neutrino 
problem should be treated under the MSW theory. 

Because of the jump at the resonance point, the — state, emitted in region I and moving 
towards the Earth, becomes a mixture of the — and + solutions in region II (after passing 
the resonance), and of course, the same statement is true for the + solution. This jump 
introduces a complication into neutrino propagation that we describe by the concept of 
different neutrino paths. It is convenient to introduce a new two-dimensional row vector 
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index A to count the two jump possibilities. The index has two possible values: (1,-1) and 
(1, 1) corresponding to "jump" and "no jump" respectively. Then for a neutrino emitted 
into the — state, the formal product of the — and A indices gives two possible values, (— , +) 
and (—,—), which count two possible neutrino paths for this neutrino. For example, the 
(— , +) path means that a neutrino originally in the — state (in region I) is converted into 
the + state at the resonance (and stays + in region II). The same simple rule is true for 
the product of the + and A indices that counts two possible paths, (+, — ) and (+, +), for 
a neutrino that is emitted into the + state. As a result, every neutrino path propagating 
from the sun to the Earth is defined by its energy S, by its "sign" index ± that indicates the 
neutrino state at the position of emission (and in region I), and by the row vector index A 
that indicates whether the neutrino is converted [A = (1,-1)] or not converted [A = (1,1)] 
from one state to the other at the resonance. We also use a two-dimensional row vector 
index A that is the product of the indices ± and A and therefore has four possible values, 
(—,+), (—,—), (+,— ) and (+, +) corresponding to four possible neutrino paths (see Fig. 
We use the path index A as well as the ± and A indices. Note, that the ± index coincides 
with the first component of A, while A is the formal ratio of A to this first component (i. e. to 
the appropriate value of the ± index). 

If a neutrino is emitted into the + state (or into the — state) with amplitude one and 
the jump probability at the resonance is 5 2 , then the evolution of this neutrino on its way 
towards the Earth is represented by the following equation 

|^ ± (region I) > — ► b x |^±a 2 (region II) >, (3.6) 

A 

where the sum should be taken over two possible values of the index A. Here the subscript 
±A 2 is the product of the ± index and the second component of the vector index A (i. e. A 2 
- the second component of A). The absolute values of the amplitudes 6a are easy to derive 
from the jump probability formula (|3.5| ) 



16(1.-1)1 = 8. 



(3.7) 



It is also obvious that the following equivalence is correct 

EN 2 = L (3.8) 



A 



Equation (|3.8|) simply expresses the neutrino flow conservation law. 

As we saw in the previous section, the + and — states are orthogonal solutions of the 
Dirac equation at any point out of the spherical resonance region (i. e. in regions I and II). 
An eigenfunction has to be an exact solution valid in the resonance region as well. On 
the other hand, any such solution must be a linear combination of the + and — states 
out of the resonance region (in regions I and II). Two orthogonal eigenfunction solutions 
are those represented in Eq. (|3.6j ) [and shown by the solid and dashed lines in Fig. [|. 
It is convenient to take these solutions as our eigenfunctions, since they are pure + and 
— states at the point of neutrino emission. (We do not include the complex phases in 
the jump amplitudes b\ and b\ because we do not need them. See for references [fl9| , ^0|] .) 
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As a result, any propagating neutrino with a given energy can be represented as a linear 
combination of the two eigenfunctions (|3.6|) , or as a linear combination of four possible 
paths, A, with mathematically only two independent coefficients, all others being given by 
the jump formulas. 

The concept of different neutrino paths and the system of notations using indices A 
and A can easily be generalized for a case of multiple resonance crossing. The case when 
the resonance region is inside the solar core, so there are possibilities of double resonance 
crossing and no resonance crossing at all, is briefly discussed in Appendix A. 



IV. NEUTRINO EMISSION AND DETECTION 

As we have discussed in the introduction section, the correct treatment of nuclear col- 
lisions requires a careful analysis of neutrino emission and detection processes, and proper 
calculations of neutrino quantum phase, emission times and energies. We carry out such 
analysis and calculations in this section. We make use of the second quantization method, 
the method of stationary phase and an expansion in small neutrino masses to finally obtain 
the neutrino wavefunction and the neutrino time dependent flux at the Earth. 



A. Second quantization 

The charged-current interaction with emission or detection of the neutrinos is given by 
beta-decay elementary particle reactions 

p + + e~ ^z? n + u e , 

p + — ► n + e + + v e . (4.1) 

It is preferable to analyze reactions ( |4.1| ) in terms of eigenfunctions of the elementary parti- 
cles in the field of external forces while keeping the interaction Hamiltonian to be zero. These 
eigenfunctions are those for nucleons in a nucleus, eigenfunctions for electrons/positrons in 
the electric field of the nucleus, and (|3.6| ) for the neutrinos. 

Let qj be the amplitude of the j th neutrino oscillator. The index j runs over the two 
exact neutrino eigenfunctions for each energy and over all energy values. (There are two 
neutrino eigenfunctions ( |3.6| ) for each energy S. We choose them to coincide with the + and 
— states at the position of neutrino emission. As a result, a summation over these eigen- 



functions is equivalent to a summation over the ± index.) Under second quantization pl|j22 
the eigenstates of the j th oscillator are given by u nj (qj), where rij = 0, 1 represents occu- 
pied and unoccupied j th state. We imagine these wavefunctions normalized over a large 
one-dimensional box with a "volume" L. Thus, the second quantized eigenfunction of the 
neutrinos is 

4 > n\,...rij,... W ni " " " U n . • , (4-2) 

and it corresponds to rij neutrinos in the j th state. The general Schrodinger wavefunction is 

<■> X "n O,, (4.3) 

{ni,...rij,...} 
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Let us consider a process of electron neutrino emission by a single nucleus, and of electron 
(or of positron) emission if any [see reactions (4.1)1. Initially all neutrino states are empty 



a 0i, 
a n 1: 



1. 





if Enj>0. 

j 



(4.4) 



In our calculations we neglect relativistic effects. If the eigenfunctions of the proton, the 
neutron and the electron are ip p , i/j n and ip e respectively, the interaction Hamiltonian operator 
'Hint for the creation of an electron neutrino is 



n mt = G F J %$ e **J* Ve dL, 



(4.5) 



where the eigenfunctions of Fermi particles are considered as scalar field operators and G F is 
the Fermi constant of interaction. The neutrino wavelength is much larger than the nuclear 
size, so we use the contact interaction approximation. Then the matrix element for the 
creation of a neutrino into the j th state is 



/^0;....l ; .... 



n 



int 



Oa : ....a,.... $i) = G F (il)j\v e ) x / if> p ip e ijj*dL = (ipj(X N ,t)\u, 



(4.6) 



where and are the second quantized wavefunctions for the initial and final states of 
the proton, the neutron and the electron, ipj is given by Eqs. (|3^ ), (|2.13|) and (|2.14|) , X N (t) 
is the position of the nucleus, and we denote the product of all factors, which depend on the 
neutrino energy and are the same for the + and — states, as j£. 

For an amplitude ao 1) .„i. > ..., corresponding to the creation of a neutrino into the j th 
neutrino state, we get 



d 

dl"»- A - 



'j '•'() ....!,:. 



Hi, 



int 

e 



(Po u ...o jr 

■iSfi t 



■Ojy 



-i£fi t 



(4.7) 



Here £fi = £{ — Sj is the change of the total energy of the nucleons and the electron in 
the frame co-moving with the nucleus during the neutrino emission, £ is the energy of the 
neutrino, and we use initial conditions (|4.4D . 

For the emitted neutrino we have at the position of emission 



\lpv>~- 



E 

± 



ue "u,.... \ipj> d£, 



(4.8) 



where 



L dk 
2^ d£ 



(4.9) 



is the density of neutrino states, and we use integration over the energy and summation over 
the ± index instead of the summation over j (remember that our eigenfunctions coincide 
with the + and — states at the position of neutrino emission). Integration of Eq. Q4.71) and 



substitution of the result into Eq. ( f4.8|) give us the amplitude of the emitted neutrino 



12 



\lpv{x,t) > 



± 



7£ c^(^(X N (t'),t')k) e - "**' |^fo*)> dt'd£. 



(4.10) 



Here X-^{t') and £' are the position and the time of neutrino emission. The density of states ui£ 
is approximately the same for the + and — states neglecting terms of the order of A/£ <C 1. 



Now we use Eqs. (|2~T|), ( |2.13| ) and ( |3 .6| ) for (ipj(X^,1f)\u e ) and for the eigenfunctions > 



to finally get the amplitude of the emitted neutrino at position Xe and time t at the Earth 



|Vv(Earth) >= -i^^] // u(S, X N ) cos6» + v (£, X N ) 

± A 

where the exponent $a is given by 

$ A = / k\ dx + £(t' -t)- £ f i t'. 

Jx n {t>) 



sin 9 

















e l * A dt'd£, (4.11) 



(4.12) 



Here index A, as usually, is the product of ± and A indices. Note that in Eq. ( [4.1 1[ ) u{£ , X N ) 
and v(£, Xj$) are calculated at the position of neutrino emission using Eqs. ( |2.13|) and ( p,14|) . 
They depend only on the index ±, i. e. only on the first component of the neutrino path 
index A. On the other hand, ue and ve are calculated at the Earth, therefore they depend 
only on the second (the last) component of the path index A and they are given by eq. ( [2.171 ). 
The exponent $a is different for different neutrino paths. It depends on both components 
of A. This is because k(x) can be either k + (x) or k-(x) for different space intervals of the 
integration in Eq. ([4.12|) . 



B. Method of stationary phase 

$A is a rapidly varying function of the neutrino emission time if and the neutrino en- 
ergy £ , so we may evaluate the double integral in Eq. ([4.1 1|) by the method of stationary 



phase [fL6 l. That is, for each of the four neutrino paths, A, we determine t' and £ by equations 

<9$a f x v dk A , . , 

A dx + (t'-t), (4.13) 



d£ Jx N (t') d£ 
0=^ = - iA «|£> +(*-*). (4.14) 

Each of the four paths, A, gives a contribution to the neutrino wavefunction at the Earth at 
detection time t. These contributions come from four different values of t' and £ (one pair 
for each path). In other words, each neutrino path has different emission time and energy. 

Equation (|4.13|) essentially yields causality. That is, since dk\/d£ is the reciprocal of the 
velocity, this equation expresses the fact that the difference between the time of emission 
t' and that of reception t should be the time of flight for the neutrino path between the 
instantaneous position of the nucleus Xj$(t') at the time of emission and the position Xe 
where the neutrino is received. Equation ( [4. 14 ) simply represents the Doppler shift of the 



energy of the emitted neutrino path due to the motion of the nucleus. Since k\(x) is different 
for neutrino paths (marked by different index A, or alternatively by indices ± and A), the 
times of emission if are different and therefore the nucleus is at a different place X^(t') and 
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has a different velocity at these times. We will solve Eqs. (|4.13|) and (|4.14|) for the four 
solutions for t' and £ as expansions in the neutrino masses. 

To zero order we neglect neutrino masses with respect to the energy. Then all four 
neutrino paths are emitted at the same time and place, to an d Xq, and with the same 
energy £q (in the laboratory frame), which are given by Eqs. (|4.13|) and (|4.14|) to zero order 



Xe — Xn — t — to 



(4.15) 



So 



1-K 



(4.16) 



We use these solutions ( 4.15 ) and ( 4.16|) in the zero neutrino mass approximation to calculate 
all factors in Eq. ( [4.11 ) except the important exponential term exp[i<£>\(8,tf)], for which we 
still need to solve Eqs. ([4.13 ) and ( 4.14j) with more accuracy. 

An application of the method of stationary phase to the double integral in Eq. (|4.11[) 



gives us 



I W >= r(£ ) EE j a [u{So, x ) cos e + v(s Q , x Q ) sin e 



± A 



(4.17) 



where r(£o) includes all factors that are essentially the same for all neutrino paths (i.e. in- 
dependent of indices ± and A), and <3>a is given by Eq. Q4.12|) with t' and £ obtained by 
solving Eqs. ( 4.13 ) and ( 4.14j) for each path. 

The probability of neutrino detection is proportional to the square of a interaction matrix 
element that is similar to the matrix element (|4.6|) . The probability of detection of a neutrino 
emitted by one nucleus is proportional to | <u e \ip u > | 2 , where \ip u > is given by eq. ( |4.17| ). 
As a result, the neutrino spectrum, V(£q), that we detect, is 



V(£q) = Vst(£ ) {V M sw(£o 



where 



E E ^ A \u{£q, Xq) cos 9 + v(£q, Xq) sin 6 ue cos 6 + ve sin 6 

± A 



(4.18) 



(4.19) 



is the factor due to the MSW effect and Vst(£o) is the neutrino energy spectrum that we 
would detect if there were no MSW effect {0 = 0, mi = m 2 = 0, Vmsw — !)• The brackets 
(...) in Eq. ( f4.18[ ) mean that the MSW factor "Pmsw should be averaged over all nuclei, 
i.e. over their positions and their motions. In the next sections we solve Eqs. ( f4.13|) , ([4.14|) 
and we average Pmsw to obtain the detected neutrino spectrum including the MSW effect. 
Let us rewrite the MSW factor ([4.19|) using the path index A as 



V- 



MSW 



A 



EI^aI' + E E \BaM cos($ A2 -$ Al ), (4.20) 

A Ai A 2 ^Ai 



where 



u(£ , Xq) cos 9 + v(£ , Xq) sin 9 w E cos 9 + ve sin 9 



(4.21) 
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Note that b\ and B\ are complex numbers. However, we need only the absolute values of Ba 
as long as we are not interested in phases which are constant in time (see Sec. [V|). That is 
why we drop the complex phases Arg(BA 2 B Ai ) in the cosine in formula (|4.20|) . [We sacrifice 
mathematical completeness for simplicity without affecting any final results.] 
From Eq. ( |2.17| ) we obtain 

«e,+ cos 9 + t>E,+ sin 9 = cos 9, «e,- cos 9 + v e _ sin 9 = sin 9. (4.22) 

It can also be verified that 



u + (£o,X ) cos 6* + v+(£q,Xq) sin6> = cos^m, 
u_(£ ,X ) cos 9 + V-(£q,X ) sin 6* = sin^M, 

where the matter mixing angle 9m(£o,X q ) is defined by 



(4.23) 



^2,*=*!™*. (4.24) 
A cos 29 — g 

9 M is equal to 7r/2 at a density far above the resonance density (g ^> A), to 7r/4 at the 
resonance point (g = A cos 29) and to the vacuum mixing angle 9 in the vacuum (g = 0) [Q . 
Here as before A = (m| — itl\)/2£q. 

Let us calculate the MSW factor averaged over both the position of emission and the 
position of detection (i.e. over the annual motion of the Earth). In this case all cross 
terms (i.e. the terms that contain phase differences) are averaged out and we have only four 
terms left in Eq. ( |4.20|) , corresponding to the four possible neutrino paths shown in Fig. |2|. 
According to Eqs. (g/7|), (|4.21|) , Q4.22|) and (|4.23|) the absolute values of amplitudes of these 
four terms are 

|£>(+,+) | = (1 — c5 2 ) *■ ^ ■* cos 9 M cos 9, 

\ B (+-)\ = £ cos M sin 0, /aok\ 
( 9 \(i/2) l 4 - /£) J 
_)| = ( 1 — <rj sin^Msin^, 

= c^sin^MCOs^, 



and the MSW factor reduces to 



(7\isw)x ,t = (l fi (+,+)| 2 + \ B (-,-)\ 2 + \ B (+,-)\ 2 + \ B (-,+)\ 2 ) 



1 

2 + 



X 

cos20(cos20 M )x o - ( 4 - 26 ) 



This is the Parke's formula [TO]. 

Without any averaging, the right side of Eq. (|4.20| ) consists of the sum of 10 different 
nonzero terms, 6 of which are cross terms, each contains phase difference factor cos($a 2 — 
$aJ. These cross terms can potentially result in interference effects. 



C. Solutions of the causality and Doppler shift equations 



Let us continue the expansion of Eqs. ( |4.13| ) and (|4.14j) as power series in m\/£ 2 <^ 1, 
m\/ £ 2 1 and g/£ 1, while the nuclear velocity is considered finite. We need 
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f = t + t(i) + iW + ..., 
X N = X + X« + X® + . . . , 

£=£ 0+ £(l)+£(2) + .... 
For the motion of the emitting nucleus we assume 

xV> = V x tV>, X^=V x t^+a x [t^f/2, 



(4.27) 
(4.28) 
(4.29) 



(4.30) 



where V x and are the velocity and acceleration of the nucleus in the x direction at 
time to (remember that the x direction points towards the Earth). We have taken the 
acceleration to be constant, which is valid in all relative cases (see the next section). Now 
we substitute (|4.27| ) and ( |4.28[) into the Doppler shift Eq. (|4.14Q and the causality Eq. (|4.13|) 
and use formula ( |2.10| ) for k. To zero order, solutions of the causality equation and the 
Doppler shift equation are Eqs. ( 4.1 5|) and (|4.16 ) as before, while the first order solutions of 
the causality and Doppler shift equations allow us to obtain and £^ 



x E 



1 - V x JXo 
1 



m\ + m\ dk'(£ ) 



4£$ 



l-V, 



d£ 

mi + m? 





rX E 


dx R 





m\ + m\ dk'(S ) 



£oa x t^ + ^^±^ + k'(£ ,X )y 3 



4£$ 



£ a x t {1) , 



d£ 



dx, (4.31) 
(4.32) 



where the last terms are the simplified variants of these two equations under the conditions 
V x < 1, V x <C a x X E . 

We now evaluate $a as an expansion in the squared neutrino masses to the second order 
accuracy. First, we use the Taylor expansion for k(S) in £ about £ Q to write Eq. ( |4.12| ), 
accurate to second order, as 



X N (tf) 

x E 



k(£ ) + (£-£ 



dk{£) 
d£ 



dx + £{t' -t) -£ fi t' 



f^E 

/ k(S ) dx-(S- S )(t' -t) + £(t' -t)- £ fl t' 
Jx N (t>) 

fX E 

/ k(S ) dx + {£ - £ fi ) (if - 1) - £ fi t. 



(4.33) 



Note that the term that contains the second derivative of k with respect to the energy is of 
third order. To obtain the second expression for $a above we use Eq. ( 4.13|) . Now using 
expansions ( [4.271 ) an d ( 4.28|) , the zero order solutions ( |4.15| ) and (|4.16|) , formula (|2.10|) and 
Eqs. (|4730| ), we find 



-Si 



Sfito 



X E 
2 



m\ + m\ 



4£n 



+ k'(S ,x) 



dx 



o 



m\ + m 2 



—Sfito + 



X 



4^o 

m\ + m\ 

4£n 



k'(S ,x) 



dx — £q 



(4.34) 
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where we drop the last term in the first equation because it is of second order and the nuclear 
velocity is small, V x <C 1. We also neglect the (1 — factor in Eq. ( 4.31 ) for t^ l \ [Note, 

that we have also neglected relativistic effects in calculating the process of neutrino emission. 
These effects give relativistic corrections in Eq. ( f4.7|) . They also give correction terms to 
the phase (4.12). These corrections are smaller than the acceleration term a.^ 1 )] jl by the 
factor V.] 

In calculating the MSW factor ( f4.20| ) we are interested only in differences of phases 
between neutrino paths. Using Eq. ( 4.34 ) we write the phase difference between the two 
paths Ax and A 2 as 



$A 2 - $A, 



)(^0j £q) + 0a(^O; Sq)<2 x , 



(4.35) 



where 



A E 



X 

—Salt 



k' A2 (S ,x) - k' Al (S ,x) 



o <-a 2 



Ai 



A 2 



t 



Aj 



(4.36) 
(4.37) 



From Eq. Q4.31| ) we have 



t 



(i) + (i) 



A, 



Ai 



f (1) , f (1) 
r A 2 + r Ai 



A E 

X 
X E 

X 



dk' A2 (S ,x) dk' Al (S Q ,x) 



dSa 



dSn 



dx 



m l + m i + dk A 2 (So, x) t dk' Al {Sq, x) 



2SI 



dS n 



dS 



dx. 



(4.38) 
(4.39) 



0o is well-known large phase sensitively dependent on the energy, which can lead to the av- 
eraging out of the neutrino flux time variations if the detected energy band is not sufficiently 
small. (paO'x is stochastic term due to nuclear collisions, and it can totally average out the flux 
variations independently of everything else. Note, that there are two effects of the nuclear 
acceleration, which produce the additional phase difference <j) a (Xo)a x in Eq. flQgp . First, 
the additional distance X*- 2 -* that the emitting nucleus moves because of its instantaneous 
acceleration, and second, the Doppler shift of neutrino energy £ W because the nucleus has 
different velocities at different emission times. These two effects give contributions, which 
differ by a factor of two. 

Let us consider the case when the indices Ai and A 2 have different second (last) com- 
ponents, i.e. the case when in the region II, after passing the resonance, the A x neutrino 
path is different from the A 2 neutrino path. Then the main contribution to the integrals in 
Eqs. ( [4.38 ) and ( |4.39| ) is given by the interval of integration in the vacuum (even for large 
values of g inside the sun, g 3> A). As a result, using Eq. (|2.15|) for k', we have 



(1) (1) _ m l 2 ±m{ 

l r A 2 =•= r Ai I — ^2 



«I E Aw 3.8 x 10^ 



/MeV\ (m\ - m\ 



8£ 4 



(4.40) 
(4.41) 
(4.42) 
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where we drop Xq, which is much smaller than X-&. For the ranges of energy and neutrino 
masses that are interesting for us, equations (|4.40|) and ( 4.42|) are sufficiently accurate. 
However, the phase 0o given by eq. ( |4.41| ) can be very large compared to one (and it is 
independent of accelerations and emission positions). Therefore, we need to calculate the 
correction to <fto due to the integration interval in the sun in eq. ( f4.36| ). 

Equations (|420| ), (fOlp and (^35|) -(^42| ) give the solution for the MSW factor for any 
given detection time and single emitting nucleus. However, to find the actual detected 
neutrino flux we still have to average the energy spectrum (|4.18| ) over accelerations and 



positions of the emitting nuclei and to integrate the spectrum over the effective detected 
energy band. 



V. AVERAGING THE TIME DEPENDENT NEUTRINO FLUX: QUALITATIVE 

DESCRIPTION 



Before we present the exact calculations let us consider qualitatively all important phys- 
ical effects that can average out cross terms and kill any observable interference time vari- 
ations of the neutrino flux. There are three effects that can potentially wash out cross 
terms: the averaging over the effective detected energy band, the averaging over nuclear 
accelerations (over collisions) and the averaging over the neutrino emission region inside the 
sun. 

Let us consider the interference between two different Ai and A2 neutrino paths. It is 
convenient to introduce the following two energy parameters. Let 



def t 

a t = n 



t 



(i) + (i) 



he 



281 



A, 



9 

mf 



def 



-.(1)1 



X E ml 

an X e m\ — m\ 



2.6 x 10" 3 KeV 



£0 



1EI 



0.40 KeV 




(5.1) 
(5.2) 



Here an is the typical (Holtsmark normal |14[|) acceleration of nuclei defined in Eq. ( B.7|) . 
Note that all numerical estimates in this section are done for the case when the Ai and 
A 2 paths are different in the region II (and consequently in the vacuum), so we can use 
Eqs. ( 4.40 ) and ( 4.42|) . [In the next section we will see that only in the case when the 
Ai and A 2 paths are the same inside the sun and are different outside does the cross term 
between these two paths survive after the averaging over the region of the neutrino emission.] 
The physical meaning of the parameters Ot and a a is explained by the equations 



ft 



d£ 



(1) 
A 2 



f Ai I So a n (t Aa ) - (t Al ) 



a H . 



(5.3) 

(5.4) 
(5.5) 



To obtain Eq. (^3) we differentiate formula fl4.36|) with respect to the energy So and use 
Eq. ( [4.381 ). Equation ( [5.4D follows directly from formula ( |4.32j ). 

First, note that the phase difference 0o varies significantly with 8 over the energy interval 
a t , so, if we integrate the neutrino spectrum over an energy band that is larger than a t , the 
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cross term between the Ai and A 2 paths is averaged out. Second, if a a /a t is large, then 



a | an is large and the phase difference ( 4.35 ) has a large random component due to the 



distribution of nuclear accelerations, a, in magnitude and direction. In this case the cross 
term is also averaged out. 

Let our effective energy band width for neutrino detection be a. It can be either the 
energy width of the spectral response function of our detection device (as in the case of the 
continuous neutrino spectrum) or it can be the thermal width of a neutrino line. We thus 
have the following two qualitative criteria for the cross term not to be averaged out 

*<« or ,<2.6xlO-»Kev(^) (^-^), (5.6) 

/MeV\ 3 (ml -m 2 \ 2 ^ , . 

a„<a t or 150 (—) (^) < 1. (5.7) 

Here we use formulas ( |5.1| ) and ( |5.2| ) for numerical estimates. The second criterion has 
another simple physical explanation: using Eq. ( |5.5|) , we rewrite the condition a a < a t as 

a H 



[i 1 ]] 2 - £ A, (5.8) 



where X = 1/8 is the neutrino De Broglie wavelength. We see that the additional distance 
that the emitting nucleus moves because of its instantaneous acceleration must be less than 
the neutrino De Broglie wavelength for the cross term not to be averaged out. However, the 
collisional decoherence is also produced by the Doppler shift of the neutrino energy (neutrino 
paths have different energies because the emitting nucleus has different velocities at different 
emission times). 

Figure |3] illustrates criteria Q5.6p , ( 577] ) in a convenient form. This figure is plotted for 



the relevant case, when the A x and A 2 neutrino paths are different in the vacuum (see the 
discussion in the next subsection). The thick solid line in the figure corresponds to the 
equality a a = a t . The criterion (|5.7|) is fulfilled in the region of parameter space below this 
line. In other words, above this line collisions destroy any coherence between neutrino paths. 
The dotted lines correspond to different constant values of at. According to criterion ( |5.6| ). 
for every given (or chosen) value of the effective detected energy band width a we must be in 
the region below the corresponding dotted line at = a for the cross term not to be averaged 
out with the integration of the neutrino spectrum over this energy band. For example, for the 
0.862 MeV 7 Be solar neutrino line, a is actually the line Doppler broadening energy width, 
which is about 1 KeV, and we find from Fig. that we can observe significant neutrino flux 
variations only if m 2 , — m\ < 3 x 10~ 9 eV 2 . 

We see from Fig. [I] that unless the effective detected energy band width a is smaller 
than ~ 0.1 KeV, the corresponding dotted line at = <J lies below the thick solid line a a = a t , 
and therefore collisional coherence criterion ( |5.7|) is redundant, and only criterion (|5.6|) is 
necessary. Alternatively, above the thick solid line the energy averaging criterion a < at is 
redundant. In particular, this is true for the detection of the 0.862 MeV 7 Be line. 

The vertical thin solid line shows the range of values of m\ — m\ that correspond to 
the neutrino flux variations of 10% and more of the standard flux (optimized over 9) in 
the case of the 0.862 MeV 7 Be solar neutrino line (we call standard the flux that we would 
observe if the standard theory would be correct, i.e. if m\ = m 2 = 0). The slanted thin solid 
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line corresponds to 10% neutrino flux variations (normalized to the standard flux) for the 
continuous solar neutrino spectrum and for the optimal choice of the detector energy band 
width and the most favorable value of 9 (see Sec. |V1B| ). Below this line the variations could 



be larger than 10%. Both thin solid lines are based on the exact results derived in Sec. [V] 
(for the dependence of these results on the mixing angle see Figs. |]-|8| and the discussion 
below). Note that the 7 Be vertical line and most of the slanted line lie below the thick solid 
line, a a — at, that is in the region where collisional decoherence effects are not significant. 

The dotted-dashed line in Fig. ^ shows the maximum possible value of m| — m\ for 
the resonance to exist inside the sun. This line is drawn for the extreme case when the 
resonance electron density n e res [see Eq. (|5.9| )1 is equal to the central solar electron density 
and the neutrino mixing angle is small. We see that for the region of parameter space 
where the collisional coherence criterion ( |5.7| ) is satisfied (below the thick solid line), the 
electron density inside the neutrino emission region (inside the solar core) is much larger 
than that at the resonance radius, and neutrinos cross the resonance region only once. If we 
suppose collisional coherence to exist, then condition (|5.7|) gives us the upper estimate for 
the resonance electron density 

(m| — m?) cos 29 



/MeV\ 

■—^ ^6.6iV A cos2# — — 




< 0.56iV A cos2fl fjj^r) cm" 3 , (5.9) 

while the central electron density is about IOOA/a- 

Finally, the dashed line in Fig. |^ corresponds to the condition Vr\t + — £_ | = Vr|t+' ) — t- \ = 
d, the inter ion distance. Here Vr is the thermal velocity of nuclei. Our approximation of 
constant nuclear acceleration is valid in the region of parameters that is below the dashed 
line. In our calculations we have considered the emitting nucleus to be accelerated uniformly 
during the time difference \t + — t-\ with the instantaneous acceleration produced by the 
nearest ion and described by the Holtsmark distribution. In fact, this model of constant 
nuclear acceleration allows us to successfully apply the method of stationary phase and to 
neglect high order derivatives in the equations of nuclear motion. We see that in predicted 
region of collisional coherence our calculations and results are consistent with the assumption 
of constant nuclear acceleration. This assumption involves only incomplete nuclear collisions, 
which are generally sufficient to decorrelate the neutrino paths. 

In Sec. ^T| we will see how criteria Q5.6Q , ( |5.7|) are related to the exact analytical calcu- 
lations and numerical results for the detected neutrino flux. However, we first discuss the 
important averaging of the neutrino flux over the neutrino emission region inside the sun. 



A. Averaging over the position of neutrino emission 

In this subsection we average the neutrino emission over all the positions of the emitting 
nuclei. In carrying out this average it is possible to consider the more complicated case 
where the resonance region is inside the solar core and the neutrinos can cross the resonance 
zero, once, or twice, as well as the simpler case where the resonance is outside the core and 
all neutrinos cross it once. 
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The emission of all neutrinos is concentrated in a region of the solar core, which has 
the radius -R CO re ~ O.2i? . Because of the high electron density in the core the energy of 
interaction between neutrinos and electrons is strong, and the following estimate is valid: 

g(0)R COTe > g(R coie )Rcare « 1500 > 1. (5.10) 

Let us consider the interference of the Ai and A2 neutrino paths. Let these paths be different 
in the region of the solar core, i.e. one of them is the + state and the other is the — state 
or vice versa. Then, according to Eq. (|2.12|) , the difference of their wave numbers inside the 
core is of the order of the larger of g and m 2 (so even for very small neutrino masses it is 
significant). Now note that the quantum mechanical phase difference $a 2 ~~ ^Ai between 
the Ai and A2 neutrino paths is given by Eq. ( [4.35D with the phase 0o defined by Eq. ( [4.361) . 



The estimate ( [5 . 1 0| ) shows that the phase O an d, hence, the phase $a 2 — $a 1 vary strongly 
with the position of neutrino emission X . In this case, the cross term between the Ax 



and A 2 neutrino paths given by formula ( [4.20| ) is the product of the smooth function of 
Xo, 'Pst\Ba 1 B/^ 2 \, and the fast oscillating function of X Q , cos[$a 2 — "^aJ- The Riemann- 



Lebesgue lemma ||16[ guarantees that this cross term is averaged out with the integration 
over the position of neutrino emission. (We give the mathematical proof of this statement 
for the case of a single resonance crossing in Appendix B). 

As a result, the cross term between the A x and A 2 neutrino paths is averaged out unless 
these paths coincide everywhere in the region of high electron density. Therefore, two 
important conclusions are valid: First, all cross terms are averaged out in the case of a 
double resonance crossing, when the spherical resonance region is in the region of the dense 
solar core. Because even if the paths only differ beyond the second resonance crossing, a 
significant portion of the last part of the paths is still in the high electron density region. 
For the same reason, the single cross term is also averaged out for the case of no resonance 
crossing at all (when neutrinos are emitted outside the resonance radius). Second, for the 
case of a single resonance crossing, there can be only two cross terms that are not averaged 
out over the region of neutrino emission. They are between the A = (+, +) and A = (+, — ) 
paths, and between the A = (— , — ) and A = (— , +) paths. 

Now note, that the cross term between the (+, — ) and (+, +) neutrino paths is 
much smaller than that between the (— ,+) and (— , — ) paths. Indeed, according to 
Eqs. Q4.20| ), ( |4.25| ) the first cross term contains the factor cos 2 #m, while the other has the 
factor sin 2 9m, and all other factors are the same. The collisional coherence criterion can be 
satisfied only if the electron density inside the neutrino emission region is much larger than 
that at the resonance radius [see the previous section]. In this case g ^> A and the matter 



mixing angle 8m is very close to tt/2, so that cos#m ~ and sin#M ~ 1 [see Eq. ( |4.24| )1. In 



other words, neutrinos are emitted mainly in the — state, which is very close to the electron 
flavor in the region of high electron density [see Fig. [I] and also Eqs. (p.21|) ]. As a result, 



only one cross term is potentially important for the neutrino flux variations, namely the 
cross term between the (— , +) and (— , — ) neutrino paths. Henceforth, we again consider 
the case when the solar core is inside the region I and we include only this cross term in 
Eq. (|4.20|) for the MSW factor. We can use Eqs. ( |4.40| ) and Q4.42|) and all other numerical 
results obtained in the previous section for this cross term, because the (— , +) and (— , — ) 
paths are different in the vacuum. 
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VI. AVERAGING THE TIME DEPENDENT NEUTRINO FLUX: 
QUANTITATIVE RESULTS 



In the next two sections we derive the exact analytical formulas and numerical results 
for the detected neutrino flux. We consider two separate cases: that of a neutrino line 
and that of a continuous neutrino spectrum. We have to treat these two cases separately 
for two reasons. First, the acceleration of the emitting nucleus leads to slightly different 
energies for different neutrino paths [see Eq. ( |4.32j| . Our detector may have different energy 
response at these different energy values. Second, the characteristic energy band width a, 
of detection, can have different origins. In the first case of a neutrino line, a is equal to 
the line thermal broadening energy width ~ 1 KeV, and we consider the energy response 
function of our neutrino detector to be constant throughout the line energy profile (i.e. we 
do not resolve the line). On the other hand, for the continuous spectrum, a is a feature of 
our detector. In this case, to get a definite numerical answer, we choose detector's energy 
response function to be Gaussian with the width 2a. Of course, in both cases we can observe 
significant time variations of the detected neutrino flux only if the criteria o a < a t and a < a t 
are satisfied. (In addition, for the continuous spectrum we must also have a a < <r, else the 
different energies of different neutrino paths will not be simultaneously included in detector 
band width.) These two criteria were qualitatively obtained in Sec. |V| from a simple physical 
consideration. We now present the accurate calculations and exact results. 



A. The neutrino line spectrum 

Let us first explore the case when the detected portion of the neutrino spectrum is a 
narrow line, as for the case of 7 Be or pep neutrinos. Actually because of the high tempera- 
ture in a region of the solar core (T« 1 KeV), neutrino lines are shifted and broadened in 
energy [p3|,[|]. They have non-symmetric energy profiles with characteristic widths approxi- 
mately equal to 1 KeV. 

Let us specifically consider, as an example, the 0.862 MeV 7 Be solar neutrino line. Its 
energy profile is non-symmetric and rather complicated For the sake of simplicity 

and for final results to be comprehensive, we fit the energy profile by Gaussian function with 
the half width e = 0.67 KeV. In other words, the standard neutrino spectrum (without the 
MSW effect) is 

V S T(£)^ex V [-(£-£ c ) 2 /e 2 }. (6.1) 

Here the central energy is S c = 862.27 KeV. We consider the energy response function of 
our neutrino detector to be constant over the whole line energy profile, so formulas Q4.18 )- 



Q4.21Q are valid with the phase difference given by Eq. (|4.35 ). We first average the neutrino 



spectrum over nuclear collisions (accelerations), and then over the positions of neutrino 
emission. Finally we integrate the detected spectrum over the energy profile ( |6.1| ). 

Let us first average the MSW factor (|4.20|) over the distribution of nuclear accelerations a. 
We need to average the cosine of the phase difference between the Ai and A 2 neutrino paths 

(cos($ A2 - $ Al )) a = (cos(0 o + 4> a a x )) a , (6.2) 
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where we use Eq. ( |4.35|) , and 0o and <p a do not depend on the acceleration. Since the 
direction of the nuclear acceleration n„ has a random distribution, we have 



(cos((p a a x )) na 
(sin(0 a a x ))„ 



1 

2 Jo 

-I 

2 Jo 



cos(0 a a cos $) sin •& dd 



sm((p a a 



4> a a 



sin(0 a a cos i?) sin $di9 = 0. 



(6.3) 
(6.4) 



Here a is the absolute value of the acceleration and •§ is the angle between the acceleration 
and the x direction, so a x = acosi?. As a result, 



(cos($a 2 - $ A J) a = cos(0 o ) 



sm[<p a a 



4> a a 



(6.5) 



We still must average ( |6.5| ) over the absolute value of the acceleration, a. For this we 

(6.6) 



use the Holtsmark distribution /h of the instantaneous nuclear accelerations JT4 

h(a)=H(P)/a H , 



where H{(3) is the Holtsmark function and (3 = f a/an- Here 



a H 



V15/ m N x ' 



(6.7) 



is the normal acceleration, Qn and mN are the electric charge and the mass of the emitting 
nucleus, (^Qf 2 ^ is the average of the three halves power of the charges of the ions producing 
the acceleration, and rij is the ion density. The Holtsmark function is 



H(/3) = — / xsmx exp —(x/j3) 3 ^ 2 dx. 

Tip J 1 J 

In Appendix C we show that the average in Eq. ( |6.5p reduces to 

sin(0 a a H ^) 



(6.8) 



sm(0 a a 



4> a a 



5> a a H /3 



H(0) d/3 = exp(-|0 a a H | 3/2 ) 



(6.9) 



Combining formulas (|4.20[) , (|6.5| ) and ( |6.9|) , we obtain the cross term between the Ai and 
A2 paths averaged over nuclear collisions 



(|5 Ai 5a 2 | cos($a 2 - $Aj) a = \B Al B A2 \ cos(0 o ) exp (-|0 a a H | 3/2 ) • 



(6.10) 



We see that unless |0 a an| < 1, the cross term is averaged out over nuclear accelerations in 
accordance with the results obtained in Sec. [V] [see Eqs. (|5.5|), (|5.7|)]. From Eq. ( 4.42|) 



8^ 4 



75 



( 



MeVV 



V £0 



1 



nin — m 



12 a v 4 ' 



10- 12 eV 



(6.11) 



where a H is calculated at the center of the sun. (|0 a a H| is roughly independent of the position 
of neutrino emission.) 
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f r i f r i 

/ ^(_ i+ )(^0) a; ) — k', J£ , x) dx = / k' + (£ , x) — k'_(£ , x) dx. (6.12) 

J Xn J X tcs 



To carry out the average over the position of neutrino emission we first note that only 
one cross term, between the A 2 = (— ,+) and Ai = (— , — ) states, is important and that 
the electron density inside the neutrino emission region is much larger than that at the 

resonance radius (see previous section). In the integral in Eq. fl4.36|) &(_+) and k'^ , are 

different only in region II. Hence, we rewrite the phase 0o as 

rX-E r . -, fXv 

As a result of Eqs. (|4.20| ), ( p.lOj ), we can write 

f V M sw(£o,t)) = (Vusw(£o)) t + A(S ) cos [M^o,t)]. (6.13) 

The MSW factor averaged over time, (Pmsw)^ is given by Parke's formula, eq. (|4 . 2 6|) , with 
cos2#m = — 1 (since the solar core density is much higher than the resonance density) 

(V M sw) t = l-(l-S 2 ) cos29. (6.14) 



2 V 2 

The amplitude of variations of the MSW factor, which we denote as A, is 
A(£ ) = 2 (\B { ^ +) B M \) Xo exp (-|0 a a H | 3/2 



5(1 -5 2 ) sin20 exp -l^anr 7 ■ (6.15) 



The factor 2 in this formula arises because we count cross terms in Eq. ( f4.20| ) twice. We 
use Eqs. ( |4.25|) with — 7r /2 to obtain the final expression for A in Eq. ( |6.15| ). Note that 



it is not necessary to average the phase 0o m formula (|6.13|) over the position of neutrino 

emission. [It is true that X res , 5 and the complex phase of the factor B^^Bf -p which we 

dropped in Eq. ( |4.20|) , depend on the position of neutrino emission X because neutrinos 



emitted at different positions in the region of the solar core cross the resonance in slightly 
different places. However this dependence is negligible for the relevant range of parameters 
where collisional coherence conditions ( |5.7| ), ( |5.9| ) are satisfied. In fact, we see that there is 
very little dependence of the cross term on the position of emission.] 

Now the only thing left to do is the integration over the line energy profile. Since the 
width of the energy profile is narrow, we consider all the factors in formula ( |6.13| ) to be 
constant over the energy range except the phase 0o(^-o^)- The Taylor expansion for this 
phase about £ c is 

M£o,t) = o (£c,t) + ^(£ - £ c ) = M£ c ,t) - - t w )(£ - £ c ) 

= M£c,t)-^X E (£ -£ c ), (6.16) 

where we use Eqs. ( |4.36| ), ( |4.38| ) and fl4.40|) [the second order Taylor term is never important]. 
Then cos 0o averaged over energy profile ( |6.1| ) is 



cos 0(£o)) £o = I J ^ V ST (£ ) cos [fofa, t)) d£ j/l J ^ V ST (£ ) d£ ^ 

= exp [-{X E Ae/2£ c ) 2 ] cos [<t> Q {S c ,t)\ = exp [-{e/2a t ) 2 } cos [<t> Q {S c ,t)\, (6.17) 
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where a t is defined by ( |5.1|) with £ = £ c . We see that unless a condition e < a t is satisfied, 
i. e. the line energy width is more narrow than a t , the cross term is averaged out over the 
energy profile. Comparing this result to criterion ( p.6| ), we see that the characteristic energy 
band width a for a neutrino line is approximately equal to the Doppler broadening width 
as one expects. 

Let F and F$t be the neutrino fluxes detected with and without the MSW effect respec- 
tively. Then their ratio according to Eqs. ( |4.18D , ( |6. 13| ) and ( |6.17| ) is 



F/FsT = (Vusw{Q) t + A(£ c ) exp [-{e/2a t f} cos [0 O (&,*)], 



(6.18) 



where we drop the constant phase tan -1 (A 2 /Ai). The eccentricity and the period of 
the Earth's orbit are e = 0.0167 <C 1 and P = lyr respectively, then Xe — (Xe) ~ 
e (Xe) cos (2irt/P), where (Xe) is the mean annual distance between the Earth and the 
sun. The time dependence of (p (£ c ,t) is easy to obtain from Eq. (|6.12| ) 



<f>o(£ c ,t) 



x E (t) 

dx + I A • dx 

(X E ) 



k+(£ci x) k_(£ Cl x) 
= (0o> + e(X E )A-cos(27rf/P), 

where we use Eq. (|2.15| ) for the difference of k' + — k'__ in the vacuum, and 



(0o) 



X r . 



kA£ c ->x) k_(£ c ,x) 



dx. 



Then the ratio of the neutrino flux to the standard flux ( 6.18Q finally reduces to 

t 



F/F ST = (V M sw(£ c )) t + A MSW (£ C ) cos 



5 ) + Ccos (2tt-) 



where 



and 



/ o\ 1 / 2 

^4msw = 5(1 — 5 ) sin26 ) exp 



C = e (X E ) A w 6.3 x 10 s 



/MeV\ 

\~£VJ 



3/2' 



exp 



ml 



10- 6 eV 



> 1. 



(6.19) 



(6.20) 



(6.21) 



(6.22) 



(6.23) 



We see that provided that ( > n/2 the factor A M sw gives the actual value of neutrino 
flux variations normalized to the standard flux, which is about 4.7 x 10 9 cm _2 s~ 1 for the 
0.862 MeV 7 Be line. If C ~ !, i-e. m\-m\< 1.6 x 10~ 10 eV 2 £: Me v, then the variations are 
suppressed by a factor ( times another factor of order unity that depends on the fractional 
part of uncertain number (<fi )/2ir. In this case the phase (0o) can be treated as a parameter 
that must be matched with observations. Particularly, if £ <C 1, then F/Fst is a constant 
in time, which can differ from the Parke's result! 

How large are the flux variations? Let us assume that m\ 3> m\ and use Eq. ( |3.5| ) for 
the jump probability and Eq. (|6.11| ) for the acceleration phase factor a flH- In this case 
^4msw, the amplitude of neutrino flux time variations normalized to the standard neutrino 
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flux [see Eqs. ( |6.21| ) and ( 6.22|) 1, depends only on two parameters: m 2 = f m 2 and the vacuum 
mixing angle 9. Figure |4] represents contours Amsw — const in the m 2 and sin9 parameter 
space for the 0.862 MeV 7 Be line. These contours are based on the exact formula ( |6.22 ). 
The maximum possible value of Amsw is 50%, which occurs for sin 29 — > 1 and the jump 
probability 5 2 = 1/2. We see that there are no significant time variations of the detected 
neutrino flux unless the neutrino mass is quite small, m 2 ~ 5 x 10" 9 eV 2 , and the mixing 
angle is quite large, sin6> ~ 0.2. Figure ||(a) an d H(b) support this statement. These two 
figures show the maximum possible value of m 2 and the minimum possible value of sin^ 
versus the desired level of the relative (to the standard flux) neutrino flux variations Amsw- 
If we neglect nuclear acceleration, i. e. set an to be zero in formula ( |6.22| ), then all figures 
will not noticeably change. This is not surprising, because the acceleration averaging factor 
in formula ( |6.22| ), exp [— |0 a aH| 3y/2 ], starts to play its role only when the energy averaging 
factor, exp [— (e/2a t ) 2 }, is exponentially small. Thus, acceleration averaging is not significant 
for the neutrino line. However, the situation is completely different if we look inside the line 
(see the next subsection). 



B. The continuous neutrino spectrum 



In the case of the detection of the continuous neutrino energy spectrum the detected 
energy band is determined by the energy response function of our detection device and, as 
we know, must be quite narrow to satisfy criterion ( |5.6| ). Moreover, formulas ( t4.18| )-( [4.21| ) 
must be modified since the detector energy response function can no longer be considered 
constant over the observed energy band. Actually, the detector response is different for 
neutrino paths with different indices A because they have slightly different energies given by 
Eq. ( ET32D . As a result, we now must allow in our calculations for the energy dependence of 
our detector response. 

We consider a neutrino detector with a simple Gaussian energy response function 



S(£) = S c exp 



(£-£ c 



a 



(6.24) 



where £ c is the energy at its center, and a -C £ c is its effective energy width, which we may 
choose in order to maximize the time variations of the detected neutrino flux. (Note, that 
the variance of S(£) is cr 2 /2.) The detection rate is proportional to the square of the product 
of the electron neutrino amplitude and the square root of the detector response function, so 
that the observed neutrino flux F is 



oo 

F(£ c , a) = P S t(£c) { [E S I(£c) ] / 3{So) d£ + 



+ E E 1^(4)^(^)1 (/ y/S(£ A2 )S(S Al ) cos($ A2 - $ Al ) d£ ) ax \, (6.25) 



Ai A 2 ^Ai 



where we consider all factors to be constant throughout the detected energy band except 
the phase difference $a 2 — $Ai and the detector energy response function S(S) . [Note that 
the energy response function should be calculated at the exact value of neutrino energy 8 a, 
which is different from £ by the value (|4.32|) .1 
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Let e d = £ - £ c < £ c . Using Eqs. ([O^), flE37D , ( ^35| ), flBTTBD and ( ggg ) we have 

OO 

/ ^S(£a 2 )S(£ Ai )cos($a 2 - $aJ dfib 



OO 7T OO 



exp 



'2^ 



-oo 



(e + £ t fa a cos #) + (e + Ai a cos ^) 

x cos |0 O (4) - e[t { H - t^} - ^ [(*2) 2 - (tg) 2 ] acos^|/ H (a) da i sintf dtf de. (6.26) 

Here again d is the angle between the nuclear acceleration and the x direction, a x = acos'd, 
and /h(o) is the Holtsmark distribution of nuclear accelerations ( 6.6)— ( 6.8). Introducing 



new variables /x = cos ■d and x = e + £ [t M > + t A >/V 2 > 

we get 



dcf 



,(1) ^ ,(1), 



CO 1 OO 



I ^S(£ M )S(£ Al ) cos($ A2 - $ Al ) dS ) = y J J J exp[-p 2 a 2 af i 2 /4a 2 } exp[-x 2 /a 2 } 

— oo — oo —1 

x {cos [4>o(£ c )} cos l(tfa ] - t M ] )x] + sin [<f> (£ c )] sin [(*5S ~ t&x] } ^(a) rfa ^ dx > ( 6 - 27 ) 
where we use definition (|5.2j ) for cr a , and as before /3 == a/an. The integral over x can easily 



be done (assuming that a a is independent of x -C £ c ) to obtain 

oo 

J \/S(£a 2 )S(£ Ai )cos($a 2 - $Ai) d^o) 

— oo 

= tt 1/2 S c a exp [-<x 2 /4(T 2 ] cos [MQ\ J J ex P MV^/ 4 ^ 2 ] H(/3) d/3 dfi. (6.28) 



1 oo 







Here we use Eq. ( |6.6| ), and a t is defined by Eqs. ( [5.11 ). The integral above with H(f3) given 
by Eq. (|6.8|) is calculated in Appendix D. The result is simply 



oo 

J ^S(£ A2 )S(£ Al )cos($ A2 - $ Al )d£ Q 







oo 






3 4 " 2 


J xexp 














dx. 



(6.29) 



As before, there is only one important cross term, that is between the A = (— , +) and 
A = (— , — ) neutrino paths. Combining Eqs. ( |6.25| ) and ( |6.29p we have for the neutrino flux 

F(£ c ,a) = 7r 1 / 2 S c aPs T (£c)(nisw(£c)> t + S c A(£ c ,a)P ST (£ c ) cos [0 O (£ C , t)), (6.30) 

where (Vusw)t * s gi ven by Parke's formula Q6.14Q , the phase <p (£ c ,t) is given by Eq. ( |6.19|) , 
and 
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A(S c ,a) = 8 (£(_,+)£? 



»(-,-) 



oo 






3 4 " 2 " 


y xexp 










4«7 t <5(l -5 2 ) (1/2) sin20 Z(f,a). 



Here we use Eqs. ( fl.25| ) and introduce two new dimensionless variables 



The function a) is 



* def / 

? = 

def 



(6.31) 

(6.32) 
(6.33) 



Z(£,a) = a£ 2 exp [- 



-£74] j a- exp 
o 



— X 



(6.34) 



The physical meaning of £ and a is easy to understand from criteria (|5.6| ), (5.7). The 
effective energy width a is the feature of our detector. Hence £ is a free parameter that may 
be chosen to detect the maximum flux variations according to formula (|6.30| ). If condition 
a < a t is not fulfilled, the cross term is averaged out over the energy according to crite- 
rion ( |5.6| ). On the other hand, if we choose a to be too small, then we detect only a small 
number of neutrinos in the chosen energy band. [The reason for the factor a 2 in eq. ( |6.31| ) is 
that the fraction of accelerations of the emitting nuclei which are small enough for both neu- 
trino paths to have energies inside the a band is proportional to ex.] Thus, for a given a there 
is the optimal value of the parameter £ = a /at [£ O pt(a0] that maximizes the function Z(£, a) 
[Z max (a) = Z(£ op t,at)], i.e it maximizes the detected neutrino flux time variations in SNUs. 
Figure |6| shows this optimal value of £ (the solid line) and the corresponding maximum value 
of the function Z(£, a) (the dotted line) as functions of the parameter a = o- t /a a . 

The standard neutrino flux (without the MSW effect) is simply 

F ST (£ c ,a) = ir 1 ' 2 S c a P ST (S C ). (6.35) 



Now using Eq. ( |6.19|) and definition (|6.23| ), we obtain the ratio of the observed flux (|6.30|) 
to the standard flux for the optimal choice of £ 



F/FsT = (V M sw(£c)) t + A M sw(£ c ) cos 



/> )+Ccos(27r- 



[again we now see that we do not need the constant complex phase of the factor B 



This formula is similar to Eq. ( |6.21|) , but now /Lmsw is given by 

fa) 



(6.36) 



(-.+)- B (-,-)]- 



.4 



MSW 



An- 1 / 2 5(l-5 2 ) {1/2) sin 26 ^ max 



Copt («) 



(6.37) 



Equations ( |6.30| )-( p.36| ) represent the final result for the case of the continuous neutrino 
spectrum. In the limiting case that a — > 0, i.e. a a ^> a t , we have 



a 

Copt (a) 

^max(^) 



o, 

0.451a£ 2 exp(-£ 2 /4), 
2, 



(6.38) 



0.664 a 



0. 
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and the flux variations are small in accordance with collisional decoherence criterion fl5.7|) . 
On the other hand, when a — > oo, i.e. a a <C o~ t , we neglect nuclear accelerations and we 
have 

a — ► oo, 

Zfoa) — 0.443 £exp(-£ 2 /4), 

4ax(«) — > 0.380. 



Thus, according to formulas (|6.36 ) and (|6.37 ), the maximum amplitude of the neutrino 



flux variations, divided by the standard neutrino flux passing through the corresponding 
optimized energy band width (we optimize it to observe the maximum amplitude of the flux 
variations in SNU), is 42.8%. This occurs for sin 26 — > 1 and the jump probability 5 2 = 1/2. 
Note, that we chose £ in order to maximize the relative amplitude of the flux variations, it 
would be larger. 

Let us assume m 2 = m\ 3> m 2 and use Eqs. (|3.5|) and (|6.11|) in the same way as 



we did in the case of neutrino line detection. The amplitude of the neutrino flux time 
variations normalized to the standard neutrino flux, ^msw, is given by Eq. (|6.37|) . For the 
optimal choice of the parameter £ (i.e. a) it depends only on three parameters: S c , m? 
and sin 8. Figures 0(a), 0(b) and 0(c) represent contours for the relative flux variations, 
^4msw = const, in the m 2 and sin# parameter space for three values of energy, 0.1, 1 and 
10 Mev (solid lines). Figure 0(d) shows constant energy curves, £ = const, for the fixed 
variation level v4 MSW = 10% (solid lines). To understand the effect of nuclear acceleration, 
we plot the corresponding contours for the limiting case — > 0, as dotted lines if the 
effect is noticeable. We see that nuclear acceleration makes a real difference, especially for 
small neutrino energies and for small levels of flux variation. Figure |8](a) and |8](b) show the 
maximum possible value of m 2 and the minimum possible value of sin 9 respectively versus 
the desired level of the relative neutrino flux variations Amsw for different neutrino energies. 
All these figures are based on the exact formula fl6.37| ) and are for the optimal choice for a. 



We see again that unless neutrino masses are quite low and the mixing angle is large, we do 
not detect time variations of the neutrino flux. 



VII. CONCLUSIONS 

There are three culprits that can destroy the coherence between different neutrino paths 
and wipe out the time variations of the detected neutrino flux. These are: (1) the aver- 
aging over the region of neutrino emission, (2) the averaging over the effective detected 
energy band, and (3) the collisional decoherence. We concentrate on the third one in this 
paper. However, to properly treat the nuclear collisions, we must also include the two other 
decoherence effects in our theory. 

Introducing nuclear accelerations into the theory requires a careful treatment of the 
process of neutrino emission (based on the second quantization method) and correct cal- 
culations of different neutrino emission times and energies. In solving this problem we use 
the convenient concept of the neutrino paths, with the two neutrino eigenfunctions as mix- 
tures of them. Different neutrino paths have different quantum phases and interfere at 
the Earth. They also have slightly different emission times and slightly different energies 
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which are determined by the method of stationary phase, which leads to the causality and 
the Doppler shift equations. Thus, there is a correlation between the flux of neutrinos at 
two nearby energies. This result contradicts those of Stodolsky P5 |, who states that for 
any steady source of neutrinos there should be no such correlations. We believe that the 
difference arises because he neglects the time dependent two point correlation function for 
the neutrino emission. These correlations are due to coherent fluctuations of the nuclear 
densities. As is well known, such fluctuations lead to correlations at different energies. They 
are produced by the nuclear accelerations we consider. 

We calculate the neutrino flux time variations and the effect of nuclear collisions (best 
summarized in Fig. separately for the case of a neutrino line and for the case of the 
continuous neutrino spectrum. 

1. Solar neutrino line. We obtain the known results for this case. The averaging of the 
neutrino flux over the line Doppler broadening energy profile washes out any flux variations 
unless the neutrino masses are very small m\ — m\ < 5 x l(r 9 eV 2 and the mixing angle is 
quite large sin6> > 0.3. 

2. Continuous solar neutrino spectrum. In order to observe flux time variations in this case, 
either neutrino masses should be very small (depending on the energy, see Fig. §a), or the 
energy resolution of our detection device should be incredibly high, o pa a t with o t given by 
Eq. ( |5.1|) . See Fig. ^ for detailed contours of the relative flux variations. 

However, for both cases and for neutrino masses larger than a certain limit (m| — 
m\ > 10 _7 eV 2 ^mIv) nuclear collisions always destroy neutrino flux variations; and for 
masses less than this nuclear collisions are unimportant. 

Recently, Bahcall et al gave three estimates for MSW solutions that best agree with 
experiments J7]. Even for their low mass solution, m\—m\ = 7.9 x 10 -8 eV 2 and sin 6 = 0.63, 
no time variations of the neutrino flux from the 0.862 MeV 7 Be line would be detected. As 
for the "vacuum oscillation" solutions, m 2 . — m\ ~ 10~ 10 eV 2 , the collisional decoherence is 
negligible for them. 
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APPENDIX A: NEUTRINO PATHS AND EIGENFUNCTIONS FOR THE CASE 
WHEN THE RESONANCE IS INSIDE THE SOLAR CORE. 



Let us consider that the spherical resonance region happens to be inside the region of 
neutrino emission. In the simplest case a neutrino emitted outside the resonance radius does 
not cross the resonance at all. In this trivial case, we do not need to use indices A and A. 
There are two neutrino paths and two neutrino eigenfunctions which coincide with + and 
— neutrino states, and they are marked by the single index ±. 

In the more complicated case a neutrino, which has been emitted at the back side of the 
sun, crosses the resonance region twice. Because of geometrical symmetry, probabilities at 
the first and second resonance crossings for the neutrino to jump from one state to another 
are the same and equal to 5 2 . Indices A and A are now three-dimensional. The four possible 
values for A are (1,1,1), (1,1,-1), (1,-1,-1) and (1,-1,1). The index A is again the 
product of indices ± and A. So, there are eight neutrino paths and correspondingly eight 
values of the path index A: (+, +, +), (+, +, -), (+, -, -), (+, -, +), (-, -, -), (-, -, +), 
(— , +, +) and (— , +, — ). For example, the (+, — , — ) path means that a neutrino originally 
in the + state is converted into the — state at the first resonance crossing and stays in the — 
state after passing the second resonance crossing. The two neutrino eigenfunctions coincide 
with + and — neutrino states at the point of emission and they (the eigenfunctions) are 
given by the formula 

[^(position of emission) > — > b\\ip±\ 2 (between the two resonances) > 

A 

— ► 51 W±a 3 (vacuum) >, (Al) 

A 

where the sum should be taken over all four possible values of the index A. Here the subscripts 
±A2 and ±A3 are the product of the ± index with the second and the third components of 
the vector index A respectively. The absolute values of the amplitudes b\ and b\ are 



s2 



1/2 



and 



|6(i,i,i)l = i 1 

|6(i,i,-i)l = °> (A2) 

|6(i,-i,-i)l = ^ 

|6(i,-i,i)l = o, 



|6(i,i,i)l = 

f 2 > ' ' 



16(1,1,-1)1 = J, (A3) 



|6(i,-i,-i)l = S(l -S 2 ) 
|6(i,-i,i)| = 5 2 . 

Of course, the choice of 6a for the unused indices is arbitrary. 
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APPENDIX B: AVERAGING OVER THE POSITION OF NEUTRINO EMISSION 



Let us consider the case of a single resonance crossing. In this appendix we show that 
the cross term between two neutrino paths Ai and A 2 is averaged out over the region of 
neutrino emission if these paths are different in the region of the solar core. Let for example 
Ai = (— , — ) and A2 = (+, +). Since the phase 4> a a x is much smaller than the phase 0o m 
Eq. ( 4.35Q , we can drop it (but only for the averaging over the neutrino emission region. It 
still can be much larger than one, resulting in the collisional decoherence if the solar core 
paths are the same.) Then, the cross term between A 2 and paths is |5a 1 -Ba 2 | cos(0 o ) [ see 
Eq. ( |L2(J| )]. Henceforward we use dimensionless variables and functions (we use the same 
notations for them as for dimensional variables): 



r = r/R Q 

x = x/Rq 

y = r±/R Q 

n e = n e /n e (0) 
A = A/g(0) 
P(r) 



— radius-vector normalized to the radius of the sun, 

— dimensionless x-coordinate, 

— dimensionless perpendicular radius r± = \Jr 2 — x 2 , 

— electron density normalized to that for the center of the sun, 

— A = [m\ — m\)/28 normalized to g(0) = \/2G F ri e (0), 

— normalized neutrino production rate inside the sun. 



It is obvious that P(r) is zero outside the region of the solar core and 



P(t) d 3 r = 4tt / P(r) r 2 dr = 1. 




(Bl) 



Now using Eqs. ( |4.36|) , (|2.12|) and recalling the definition g = f y/2Gpn e , we get for the 
cross term averaged over the position of neutrino emission 



B Al B A2 \ cos( 



cos 



|r|<r c 



A / n{f) dx 



d 3 r. 



(B2) 



Here 7Y orP is the dimensionless radius of the solar core and we introduce 



^ djf gi^Re ^ 27000 > 1, 
he 

f(v) d ^\B Al B A2 l 

K {r) = [n 2 - 2nA cos 29 + A 2 

Using cylindrical coordinates x and y, we write 

Y 



1/2 



-Y x 



(B3) 

(B4) 
(B5) 



dx, (B6) 



where Y 



def 



y 2 . Using a new variable t instead of x 
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def 



K(r) dx, 



(B7) 



we rewrite Eq. ( p6|) as 

(\B Al B A2 \ cos(0o)) r 

where 



Re / 2iry dy 

T. 



P(r) f(x,y) 
k(t) 



exp [— iAt] dt, 



(B8) 



T ± = - 



Kir) dx. 



±Y 



Integrating Eq. ( |B8|) by parts and using equation P(r 
outside the region of the solar core), we get 



(B9) 

(no neutrino emission 



' core 

|S Al -B A2 |cos(</)o)) r = Re / 2nydy / — - — - 

J J lA Kir) 



1 x d 



(r) r dr 



P{r)f{x,y) 



Kir) 



exp [-iAt] dt. (BIO) 



where we also use d/dt = (1/ n)(x/r)(d/ dr) . To estimate the absolute value of the averaged 
cross term we drop the exponential factor in Eq. (|B10| ) and we return to the old variable x 
instead of t. We have 



fcore Y 

\(B Al B A2 cos(0 o )) r | <- J 2nydy J 

—Y 



d P{r) f{x,y) 



dr 



Kir) 



dx. 



(Bll) 



We see that now we have a small factor 1/A ~ 4 x 10 -5 in the formula for the averaged 
cross term. We can integrate Eq. ( [B10|) by parts again and get the factor 1/A 2 , and so on. 
However, instead we content ourselves with the numerical estimate from Eq. ( |B11| ). 

Since we consider the case of a single resonance crossing, the resonance radius is outside 
the neutrino emission region, n(r core ) > A and k(t) ~ n(r) [see definition (BS)]. Let us 
further assume that f(x, y) = \B Al B A2 \ is constant (as in the case when the electron density 
inside the neutrino emission region is much larger than that at the resonance radius), and 
replace f(x,y) by its largest possible value f(x,y) = 1/4, according to Eqs. ( |4.25| ) and for 
our choice A 2 = (+, +) and Ai = (— , — ). In this case, using the spherical coordinate system 
x = r cos a and y = r sin a, we rewrite Eq. ([Bll ) as 



' corf 

[B Al B A2 cos(0 o )) r | < j J 




7T 

2A 



d P(r) 



dr K\r) 



■k/2 



dr 



o 



d_ 

dr 



P(r) 
n(r) 



sin a cos a da 

P(r) 



dr = — 
A 



o 



n(r) 



dr w 2 x 10 



-l 



(B12) 



For the final numerical estimate we use the production rate and the electron density for the 
standard solar model HQ. We see that an averaged cross term is indeed very small if the 
initial neutrino states are different. 
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APPENDIX C: DERIVATION OF EQUATION 



Substituting the Holtsmark function (6.8) into Eq. (|6.9| ) and introducing a new parameter 
3 get 

sm(0 a a) \ 2 f f sm{Ajj) 



A = f a a H , we get 



7T 







Af3 



exp 



— (x//3) 3 ^ 2 — sinx dx d/3. 

P 



(CI) 



Now we change the variable (3 to a new variable t = f A/3 to obtain 



sin (0 a q) 
(p a a 



2 /■ sin t . 

— (it / exp 

71 J t 



7T 



A 3 / 2 (x/t) 3/2 
/ 2 dt J exp [-A 3 / 2 (x/t) 



OO OO 

2 /" Sill t 



3/2 



x sin x 
d 



dx 



fsinx — x cosx) dx. 



(C2) 



We integrate the second integral in Eq. ( |C2| ) by parts to get 



sm <z>„a 



-A 3 ' 2 J J x 1/2 r 7/2 sint exp [-A 3/2 (x/t) 3/2 ] (sinx - x cosx) dx dt. (C3) 







Introducing a new variable y = x/t instead of the variable x, we obtain 

3 



sm n a 



4> a a 



71 



X 



'-A 3/2 / v 1 / 2 



y ' exp 



sin (t) sin (yt) 



■■X. 

dt-yj 



dy 



sin (t) cos (yt) 



dt 



(C4) 



Now using formulas p4 | 



sin(t) sin (yt) 



dt 



vr/2, y > 1 



and 



sin(t) cos (yt) 



vr/2, y < 1 
dt = { tt/ 4, y = 1 , (C5) 
0, y> 1 



we get 



sm 0„,a 



a a 



A 3 / 2 / y 1 / 2 exp 



-A 3/2 y 3/2 \ dy. 



(G6) 



Changing again the variable of integration to x = y 3 / 2 and substituting A = a an, we 
obtain the final result 



sm {(p a a) 
<p a a 



A 3 ' 2 



exp 



-A 3 ' 2 



x 



dx = exp 



_ A 3/2 



exp 



.3/2 



(G7) 
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APPENDIX D: DERIVATION OF EQUATION flQ9D 



Let us introduce a parameter A = f a a /a and denote the integral in Eq. QB.28|) as 1(A), 



1(A) 



def 



1 oo 


" A 2 




// exp 


T 










H((3) dp dfi 



1 oo oo 



exp 







A 2 

4 /3V 



exp 



(x/(3) 3 ^ 2 — sin x dx d/3 dfi. 



Changing the variable x for a new variable y == x/(3, we get 

_ 1 oo oo 



7T 



/(A) = — I d/j, I y exp —y 3 ^ 2 dy I (3 exp 



o o 



s 



sin (/3y) d/3. 



Using the standard formula for the integral over j3 ||24|| , we obtain 



(D2) 



1(A) 



71 



1/2 ^3 



1/ exp 



3/2 



/ exp 



A 2 ft 2 



dfi 

V 3 ' 



(D3) 



Now introducing a new variable t = f l//i 2 instead of //, we get 



7T 



1/2 ^3 



I/ 2 ex P 



3/2 



dy / exp 



Z/ 2 ' 
A 2 * 



dt 



71 



1/2 a 



exp 



.«3/ 2 _ y_ 
y A 2 



dy. (D4) 



Changing the variable of integration to x = f y 1 / 2 and substituting A = er a /cr, we obtain the 
final result 



I(aJo) 



4 a 



7rV2 an 



oo 






3 4 °" 2 " 


J x exp 


- x — 








(D5) 



36 



FIGURES 




n e 

FIG. 1. The dependence of the neutrino relativistic mass mj? c | = £ 2 — k 2 on the electron density. 
We see that the — state is almost the electron flavor for high densities and it's almost the muon 
flavor for small densities (if the mixing angle is small). The opposite statement is true for the + 
state. Nuclei emit mostly into the — state. If the conversion is adiabatic, the resulting state is the 
muon difficult to detect at the Earth. 
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Resonance 




FIG. 2. All four possible neutrino paths are shown on the upper picture for the case with the 
resonance outside the solar core. The left bottom picture shows two adiabatic paths A = (— ,— ) 
and A = (+,+). This picture can be compared with Fig. |l[ On the right bottom picture two 
nonadiabatic paths A = (— , +) and A = (+, — ) are shown. There are two neutrino eigenfunctions, 
one indicated by the solid lines and the other indicated by the dashed lines. Each of these eigen- 
functions is a liner combination of the corresponding two paths. An electron neutrino emitted at 
the sun can still stay an electron at the Earth due to a nonadiabatic jump from the — to the + 
state. 
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0.1 1 10 

E, Mev 

FIG. 3. A graphical demonstration of the dependence of the qualitative coherence criteria 
(and of some exact results) on two neutrino parameters, the energy and — mf. The thick 
solid line corresponds to the condition a a = at- Above this line nuclear collisions destroy any 
coherence, and as a result, there is no time variations of the detected neutrino flux. Below this 
line collisions are unimportant. The dotted lines correspond to different constant values of at- The 
flux time variations are averaged out over the effective detected energy band width a for the region 
of parameter space above the line at = a. The vertical and slanted thin solid lines correspond 
to neutrino flux variations equal to 10% of the standard flux (optimized over 6) in two cases 
respectively: the 0.862 MeV Be solar neutrino line and the continuous solar neutrino spectrum. 
The dotted-dashcd line shows the extreme case when the resonance electron density is equal to the 
central solar electron density and cos 26 « 1. The model of constant nuclear acceleration is valid 
for the region of parameter space below the dashed line Vr\t+ — 1-\ = d (Vr is the thermal velocity 
of nuclei and d is the mean interparticle distance in the region of the solar core) . 
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0.2 



0.4 



0.6 



sin 9 

FIG. 4. Contours ^4msw = const are shown for the 0.862 MeV Be neutrino line. These curves 
represent the amplitude of neutrino flux time variations normalized to the standard neutrino flux. 
If we neglect nuclear acceleration, the curves will not noticeably change. 




0.001 0.010 0.100 1.000 0.001 0.010 0.100 1.000 



H MSW H MSW 

FIG. 5. (a): The maximum possible value of m 2 as a function of the desired level of relative 
time variations ^4msw f° r the 0.862 MeV Be neutrino line, (b): The minimum possible value of 
sin 9 as a function of Amsw- 
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CT a/ CT t 



FIG. 6. The optimal value of parameter £ = a /at as a function of the parameter a — o"t/o"a 
is shown by the solid line. The corresponding maximum value of the function Z(£, a) is shown by 
the dotted line. 
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(a): Energy is 0.1 MeV 
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(b): Energy is 1 MeV 
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FIG. 7. (a), (b) and (c): Contours ^msw = const are shown for observation of the continuous 
neutrino spectrum at different energies (solid lines). These curves represent the amplitude of 
neutrino flux time variations normalized to the standard neutrino flux, (d): Constant energy 
curves £ = const are shown for Amsw = 10% variation level (solid lines). The dotted lines show 
how the contours would look if we neglect nuclear acceleration. 
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FIG. 8. (a): The maximum possible value of m 2 as a function of the desired level 
neutrino flux variations ^4msw is represented for different neutrino energies, (b): The 
possible value of sin# as a function of vImsw for different neutrino energies. 
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